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Abstract 

This article is devoted to analyze control properties for the heat equation with singular 
potential —fj,/\x\ 2 arising at the boundary of a smooth domain Q C R , N > 1. This problem 
was firstly studied by Vancostenoble and Zuazua [18] and then generalized by Ervedoza [8] in 
the context of interior singularity. Roughly speaking, these results showed that for any value 
of parameters fi < fi(N) := (N — 2) 2 /4, the corresponding parabolic system can be controlled 
to zero with the control distributed in any open subset of the domain. The critical value fJ,(N) 
stands for the best constant in the Hardy inequality with interior singularity. 

When considering the case of boundary singularity a better critical Hardy constant is 
obtained, namely (an '■= N 2 /4. 

In this article we extend the previous results in [18], [8], to the case of boundary singularity. 
More precisely, we show that for any \i < /in, we can lead the system to zero state using a 
distributed control in any open subset, excepting some cases where Q changes the convexity at 
origin, in which cases the control is distributed in a semi-annulus surrounding the singularity. 

We emphasize that our results cannot be obtained straightforwardly from the previous 
works [18], [8]. 
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1 Introduction 

In this article we present some new results concerning the exact controllability of the heat equation 
with singular quadratic potential —fi/\x\ 2 . 

From the mathematical viewpoint, the study of problems with inverse-square potentials is 
motivated by models which appear for instance in the context of combustion theory [15] and 
quantum mechanics [7] . 

Evolution problems with the potential — ^/\x\ 2 have been intensively studied in the recent 
decades. Among them, we remind the pionerring work by Baras and Goldstein [2] in which they 
considered the corresponding heat-like equation with the singularity localized in the interior of a 
smooth domain SI C K , N > 2 (If N = 1 they deleted the origin so that € 9f2, where dtt 
denotes the boundary of tt). They derived necessary and sufficient conditions for such systems 

BCAM Basque Center for Applied Mathematics, Mazarredo, 14 E-48009 Bilbao - Basque Country - Spain 

2 Dcpartamcnto dc Matcmaticas, Universidad Autonoma de Madrid, 28049 Madrid, Spain 



to being well-posed. More precisely, they showed the well-poscdness holds true whenever /i < 
(N — 2) 2 /4, whereas if [i > (N — 2) 2 /4 there is instantaneously blow-up for the solution. The 
critical value (N — 2) 2 /4 is the best constant in the corresponding Hardy inequality (see e.g. 
[14], [13]). Later on, the issue of singular or degenerated potentials has been also analyzed by 
the control community. The authors in [19] analyzed the control and stabilization properties 
of the corresponding wave equation. Then they showed in Vancostcnoblc and Zuazua [18] that 
the corresponding heat equation can be controlled by a distributed control which sourrounds the 
singularity. This result has been generalized in Ervedoza [8] where any geometrical constraint of 
the control region was removed. We also mention the paper [3] and references there in studying 
the control of parabolic equations degenerating at origin. Recently, the work [8] has been slightly 
improved in [17] when studying some applications to inverse problems. In all situations above the 
authors showed that the well-posedness, control and stabilization are very much related to the 
classical Hardy inequality in which the best constant is (N — 2) 2 /4. 

In this paper we consider the heat equation with the potential — /i/|x| 2 , where the singularity 
x = is located on the boundary dil of a smooth open set fl C WL N , N > 1. This work is aimed to 
extend to the case of boundary singularity the paper [8] which provides the most general control 
results in the case when the singularity is localized in the interior of the domain. 

We point out that our results cannot be deduced straightforwardly from the case of interior 
singularity and requires an independent analysis. Our main tools concern in the Carleman es- 
timates which is the classical way to prove observability properties for parabolic systems. The 
major difficulty consists in finding proper weight functions to develop efficient Carleman estimates. 
In our case, the weights in [8], [18], are not even allowed to recover the results in the range of 
parameters /i < (N — 2) 2 /4 shown in the case of interior singularity. A proper modification of the 
weights in [8], will be done here. 

Before entering into details, let us fix some ideas. 

Let N > 2 and consider a smooth bounded domain il C K N such that G <9f2, and let w C O 
be a non-empty open set. Assume also that T > is fixed. We are interested in the question of 
controllability of the following problem 

(x,t) eilx (0,T), 

(x,t) G dCl x (0,T), 
x e fi, 

where uo G L 2 (f2) and / G L 2 {Q, x (0, T)) is a function supported in the control region uj. 

The null-controllability problem reads as follows: Given any Uq € L 2 [il), find a function 
/ G L 2 (uj x (0,T)) such that the solution of (1.1) satisfies 

u(x,T)=0, xeQ. (1.2) 



d t U - Alt - i — i2"14 = /, 

\x\ 

u{x, i) = 0, 
u(x,0) = u (x), 



In order to discuss the well-posedness and null-controllability of (1.1) we need to establish the 
proper functional framework which corresponds to the problem. The crucial role of this issue is 
played by a new critical value of fj. which determines the features of system (1.1). More precisely, 
when moving the singularity from interior to the boundary, the critical Hardy constant jumps 
from (N — 2) 2 /4 to the critical value 

N 2 

A*w := — ■ (1-3) 

This is guaranteed by the new Hardy inequality with boundary singularities (see e.g. [6], [10], [9]) 
stated in a simplified form as follows. For any fj, < N 2 /4 there exists a constant 6*1(7, 17), such 
that for any < 7 < 2 the inequality 

/ f^-dx + M / fn dx < / |V U | 2 + d(7,tt) / « 2 dx, (1.4) 
Jn \ x \ ' Jn \ x \ Jn Jn 
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holds for all u £ H^(Q). 

Moreover, an improved Hardy inequality (cf. Theorem 1.2, [4]) states that for any < 7 < 2 
and any fi < fj,jy = N 2 /A there exist Ci(7, f2), (^2(7, fi) > such that 

,2 



Ci 



/ u 2 dx + 


!\ 


In 


Jn 1 



M H 2 



1 2-7 



Vul 



etc, VueiJo 1 ^). (1.5) 



This inequality will be used in the proof of the Carleman estimates. 

The main results of this paper are very much related to the geometry of fi at the origin. In 
view of that, we say that satisfies the property (P) if 



There exists a neighborhood Vo of the origin x = such that 
(P) : { x-n<0, Vi€fflnVo, 

where ft denotes the outward normal vector. 



(1.6) 



Roughly speaking, a domain Q C K 2 with the property (P) is cither concave or flat in a small 
neighborhood of the origin, as shown in the cases CI and C2 in Figure 1. 
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Case C2 
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Figure 1: The geometry of O C K 2 is drawn in red color. 
Otherwise, a domain Q C M 2 not verifying the property (P) corresponds to Figure 2 below. 



Case CI 




Case C2 




x = 






x = 










) 




71 


) 







~y^5 




y = 



Figure 2: The geometry of O C R 2 is drawn in red color. 



A domain C M. N with N > 3 satisfying the property (P), is much more complicated to 
be drawn. However, in this multidimensional case, the property (P) may be described through 
the second fundamental form corresponding to dVl in a neighborhood of origin (see e.g. [12]). 
Generally speaking, the property (P) is related to the non-convexity at origin. In particular, a 
domain strictly convex at origin does not satisfy the property (P). 

Next, we are in conditions to state the main results of this paper. 
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Theorem 1.1 (Null-Controllability). Let SI C M. N , N > 1, be a smooth domain satisfying the 
property (P) and assume fi < /ijy. Given any non-empty open set uj C SI, for any time T > 
and any initial data uq G L 2 (S1), there exists a control f e L 2 (gj x (0,T)) such that the solution 
of (1.1) satisfies (1-2). 

Theorem 1.2 (Null-Controllability). Let SI C R N , N > 1, be a smooth domain which does not 
fulfill the property (P). We also assume [i < (in and let r\ > be a fixed constant. 

Then there exists r$ = ro(Sl) > (independent of r\), < ro < r\, such that for any time 
T > and any initial data uq £ L 2 (Q), there exists a control f G L 2 (uj x (0,T)) distributed in the 
non-empty open set u> C Q, defined by 

uj := {x G SI | r < \x\ < ri}, (1.7) 

such that the solution of (1.1) satisfies (1.2). 

Following the by now classical HUM method (cf. [?]), the controllability property is equivalent 
to an observability inequality for the adjoint system 

d t w + Aw + -^rw = 0, (x,t) G O x (0,T), 

l x ' (-\ a\ 

w(x, t) = 0, (x, t) G dQ. x (0, T), V 1 - 8 ) 

, iw(;e, T) = wt{%), x g SI. 

More precisely, when < /ijv, we need to prove that there exists a constant C such that for all 
wt G i 2 (Sl), the solution of (1.8) satisfies 



Mx,0)| 2 dx < C J J \w(x,t)\* dx dt. (1.9) 

n ux(0,T) 

In order to prove (1.9), we will use a particular Carleman estimate , which is by now a classical 
technique in control theory.Indeed, the Carleman estimate we will derive later implies that for any 
solution w of (1.8), 

\w(x, t)\ 2 dx dt < C T J J \w(x,t)\ 2 dx dt. (1.10) 

s)x(if) wx(o,r) 

Let us show that (1.10) implies (1.9). Indeed, multiplying the system (1.8) by w and integrating 
in S7 we formally obtain 



1 d 
2~dt 

From (1.4) we have that 



/ w 2 (x,t) dx = / |Vw| 2 dx- fx [ dx. 
in in in F 



1 d 
2~dt 



w 2 (x,i) dx > -Ci(7,Sl) / w 2 (a;,t)dx, 
n in 



Then we get that the function t H> exp 2Cl ( 7 ' n )* \\w{-, x)\\ 2 l2 iq\ is increasing. Then we obtain 



li 2 (n) 

r 3T/4 



/ / w 2 (x,t) dx dt > I exp- 3TCl ^>/ 2 / W 2 (.t,0) dx, 
ir/4 in 2 7 n 



/T/ 

From here and (1.10) we obtain (1.9) 



(1.11) 
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Well-posedness via Hardy inequality 

Let us now define the set 



f fo Vu 2 - u N u 2 / \x 2 + Cu 2 ]dx i 
C:=\C>0s.t. inf ^0_L! 1 ,„ , J — > 1 i. 1.12 



Of course, C is non empty due to inequality (1.4). Next we define 

Co = inf C. (1.13) 

cec 

Then, for any /i < /ijv := N 2 /A we introduce the Hardy functional 

:= / |Vu| 2 cfo-M / -^cfe + Co / u 2 dx 7 (1.14) 

which is positive for any test function due to inequality (1.4) and the election of Cq. Then we 
define the corresponding Hilbcrt space as the closure of Cq°(&Y) in the norm induced by 
Observe that for any \x < /ijv the identification = Hq(Q,) holds true. Indeed, if \i < pt^ firstly 
we have 

|M|2 >h_JL) [ | V d 2 <te-^ f u 2 dx. (1.15) 

On the other hand, from the definition of Co we obtain that there is a constant C > such that 

\\u\\ 2 H ^>C f u 2 dx. (1.16) 
Jn 

Multiplying (1.16) by Coh/{C\^m) and summing to (1.15) we get that 

I MIL >o„ [ \vu\ 2 dx, 

Jn 



for some positive constant (C^ tends to zero infinity as fi tends to zero). 

Besides, in the critical case [i = fiN, H^ is slightly larger than Hq(Q,). However, using cut-off 
arguments near the singularity (see e.g. [?]) we can show that 

N|ff M(N) > CelMljji^B^ofi, VuGfl^, (1.17) 



where C e is a constant going to zero as e tends to zero and B e (0) denotes the closure of the ball 
of radius e centered at the origin. 

Let us define now the operator := —A — n/\x\ 2 +CqI together with its domain as 

D(A„) := {u 6 Hp | A„u £ L 2 (Q)}. (1.18) 

The norm of this operator is given by 

||"lb(A M ) = IMU=(n) + ||^ u IU 2 (n) V/U < hn- (1.19) 

With these definitions, by standard semigroup-theory one can show that for any fi < (J-n the 
operator (A^, D{A^)) generates an analytic semigroup in the pivot space L 2 (Q) for the equation 
(1.1). For more details we refer to Theorem II. 1, p. 3, [17], which can be adapted in the context 
of the space H^ introduced above. 
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2 Null controllability in the case fi < iV 2 /4 



First of all, to simplify the presentation, we assume that ^ ui otherwise it is straightforward 
since the control acts locally near the singularity. We also assume that fi D B\ (0) is included in 
and Bi (0) PI Q is empty. This can always be done by a scaling argument. In the sequel we also 
consider a nonempty subset wq (1(1 w whose role will be emphasized in the next paragraph. 

In what follows we are going to justify the result of Theorem 1.1. For that we will apply 
Carleman estimates using a modified Fursikov-Imanuvilov [11] weight tp which, in particular, 
satisfies 

VV>-a;>0, VxefinVo, (2.1) 

in a small neighborhood Vo of the origin. 

In view of that, in the sequel we assume that the geometry of f2 satisfies the property (P). 
Without losing the generality, by a scaling argument we may assume that the property (P) reduces 
after normalization to 

x ■ n < 0, V x G dn, \x\<l, (2.2) 



2.1 Carleman estimates. A choice of a proper weight a 

As said in the introduction, the main tool we use to address the observability inequality (1.10) is 
a Carleman estimate. 

The major problem when designing a Carleman estimate is the choice of a smooth weight 
function er, which is in general assumed to be positive, and to blow up as t goes to zero and as t 
goes to T. Hence wc arc looking for a weight function a that satisfies: 

f a(t,x) > 0, (x,t) G x (0,T), 

I (2 3) 

] lim a(t,x) = lim a(t,x) = +oo, x G f2. ' 

When shifting the singularity from interior up to the boundary the weight in Ervedoza [8] violates 
some necessary conditions to apply the Carleman estimates; the weight ip blows up at origin and 
this violates the fact that ip is constant on the boundary. 

In view of that, we propose the weight 

' l x l \ A , 



a(t,x) = e(t)(c x -\x\ 2 i;-(^) X cf>), <t> (2.4) 



where A is a positive parameter aimed at being large and ro is a positive constant aimed to be 
small enough such that it verifies 

• f 1 1 1 ■ f 1 
rn < mm < , > mm < ; ; — - — ; , 

1 100' 2 - 7 J l|I>VU-(n) + WU- ( n)' 

( C - V /(7 - X) ( C * Y /h - 1} } (25) 

VSl^li-^+SlD^I^^; >\2\v\\D4>\ L ~ m J \ { -> 

where 7 is a fixed constant such that 1 < 7 < 2. C2 stands for the constant in inequality (1.5), C\ 
is big enough meant to assure the positivity of a. Besides, ip and 9 are defined as follows. Firstly, 
we define 6 by 

with k = 1 + 2/7. Next wc introduce ip satisfying the conditions 

4' > 1, Vie O, 

i> = i, Vie on, (2.7) 
|V^|>2C n Vxen\u , 
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where the constant Co > in (2.7) is such that verifies 

\x ■ rt\ < C n \x\ 2 , yxedQ. (2.8) 

Moreover, we assume that 

x-Vijj>0, Vied, \x\ < 1. (2.9) 

Remark 2.1. Such a function ip satisfying conditions (2.7)-(2.9) exists. Indeed, let us assume tjj 
is a smooth modification of by now the classical Fursikov-Imanuvilov weight (see [11]), with the 
following property close to the origin: 

rp = d(x)>0, V x G fi, \x\ < 1, (2.10) 

where = d(x, dfl) denotes the distance to the boundary. We remark that the distance d 
yields the conditions (2.7) since it satisfies the eikonal equation |Vd(x)| = 1 a.e. In particular, the 
assumptions (2.2)-(2.10) (see e.g. [1]) ensure the validity of (2.9). 

To end up, we can choose ip := S(ip + 1) which satisfies the constraints (2.7)-(2.8) for S large 
enough. 

Motivation for the choice of the weight a. 

Roughly speaking, in order to show the Carleman estimates below, the weight a needs to fulfill 
some necessary conditions: 

CI. Vct ■ n > 0, for all x G 90. 

C2. — D 2 <j(x)(£ i , £) > 0, for all £ € M. N , x £ il, \x\ < ro, where tq is some positive constant 
independent of A. 

C3. —Act > 0, for all \x\ > ro, x ^ uq, where ro is the same constant chosen in the condition C2. 

C4. In order to get an expression in terms of the optimal Hardy inequality, the gradient term Vct 
must degenerate at origin at least as \x\. 

The weight a = C\ - \x\ 2 ip - \x\ x e x ^ fulfills the conditions CI, C2, C4. Moreover, it verifies the 
condition C3 provided 

x-Vip>0, VxeOnVo (2.11) 

for some neighborhood Vo of the origin. The normalization by ro in (2.4) is required by technical 
reasons needed precisely in the proof of Lemma 2.6. 



2.2 Main result 



We claim that 



Theorem 2.2. There exists a positive constant K\ and Xq such that for A > Ao there exists sq(X) 
such that for any s > sq we have 



s\ 2 



Ox(0,T) 



^) X 4>e- 2s °\Vw\ 2 dxdt + s J J Be 

f2x(0,T) 

?e- 2s °\x\ 2 \w\ 2 dx dt + s 3 A 4 



dx dt 



9 3 e- 2sCT (^) 3 V \ w \ 2 dxdt 
r 



S)x(0,T) 



Ox(0,T) 



< K sX 



^0 



e-^lVwr dx dt + s' 5 A 



3\4 



\x\,3X 



^e- 2sa \w\ 2 dx dt). (2.12) 



w o x(0,T) 



w o x(0,T) 



From Theorem 2.2 we can easily obtain the observability inequality (1.9) via Cacciopoli's 
inequality. The details could be reproduced step by step as in Section 2.2, p. 12, [8]. 
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2.3 Preliminaries and useful lemmas 

Now, let us assume that w is a solution of (1.8) for some initial data wt G Hq(Q), and define 

z(t,x) =e- sa( - t ' x) w(t,x), (2.13) 

which obviously satisfies 

z(T) = z(0) = in (2.14) 

due to the assumptions (2.3) on cr. The positive parameter s in (2.13) is aimed to be large. Then, 
plugging w = z cxp(cr(t, x)) in the equation (1.8), we obtain that z satisfies 

d t z + Az + -^--z + 2s\7z-\7cr + szAa + z(sdt<J + s 2 \\7cr\ 2 ) = 0, (x, t) G x (0, T), (2.15) 



with the boundary condition 

2 = o, e anx (o,t). (2.i6) 

Let us define a smooth positive radial function a(x) = ck(|x|) such that 

r , / 0, \x\ <r /2, 
t 1W R > r , 

where ro > is selected as in (2.5). 
Setting 

Sz = Az + -^z + z(sd t (T + s 2 \Va\ 2 ), 
\xy v / 

Az = d t z + 2sVz-'Va + szAa(l + aj, ( 2 ' 18 ) 

Pz = —saAaz 

One easily deduces from (2.15) that 

Sz + Az + Pz = 0, HSzH 2 + Pz|| 2 + 2 < Sz,Az >= \\Pz\\ 2 , 

where || • || denotes the L 2 (n x (0,T)) norm and < •,• > the corresponding scalar product. 
Especially, the quantity 

I =< Sz, Az > -i||s<xzAcr|| 2 (2.19) 

is non positive. 

Lemma 2.3. The following equality holds: 

I = -2s J J D 2 a(S7z, Vz) dx dt + s J J \d n z\ 2 d n a ds dt 

fix(O.T) OOx(0,T) 

-s J J \\/z\ 2 Aa a dx dt + I J J \z\ 2 A 2 a(l + otj dx dt 

fix(0,T) Qx(0,T) 

|z| 2 Va • VAer dx dt + | |z| 2 Acr Aa dx dt 

fix(0,T) Qx(0,T) 



i yy |2| 2 (s<9 2 (i + 2s 2 9 t (|Vcr| 2 )) dxdt-2s 3 Jj \z\ 2 D 2 a{Va, Va) dx dt 

Qx(0,T) fix(0,T) 

^ a|z| 2 Aa( S 2 a t a + s 3 |Va| 2 ) dx dt - L J J 



a 2 \z\ 2 \Aa\ 2 dx dt 



Qx(0,T) Ox(0,T) 



I^Acr a dx dt + 2/xs / / y\jx ■ Vcr dx dt (2.20) 



Ox(0,T) Ox(0,T) 

where d n — n ■ V and ds denotes the trace of the Lebesgue measure on dQ. 
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Here we omit the proof of Lemma 2.3 since it may be found in [8]. It is worth mentioning 
that the upcoming computations justified by integrations by parts are done formally. However, we 
notice that the final estimates make sense in our functional framework. A priori the regularity of 
the operator := —A — n/\x\ 2 +CqI is not enough to justify the integration by parts since the 
lack of regularity appears at the singular point x = 0. This issue is presented in a detailed manner 
in [5] in the context of the wave equation with singular potential localized on the boundary. 

Now, we will decompose the term / in (2.20) into several terms that we handle separately. 

Let us define Ibd the boundary term in identity (2.20): 



hd~ J J s\d n z\ 2 8 n a ds dt. (2.21) 

anx(o.T) 

Then define as the sum of the integrals linear in a which do not have any time derivative: 
h = -2s II L> 2 cr(Vz,Vz) dxdt-s // |Vz| 2 At7 a dx dt + ^ // \z\ 2 A 2 a{\ + a) dx dt 



f!x(0,T) 



nx(o,T) 
\z\ 2 Va ■ VAcr dx dt 



nx(o,T) 
\z\ 2 A<7 Aa dx dt 



Ox(0,T) 



/IS 



-Act a dx dt + 2/xs 



\4 



x ■ Vcr dx dt. (2.22) 



S)x(0,T) 



Qx(0,T) 



We then consider the integrals involving non-linear terms in a and without any time derivative, 
that is 

I ni = -2s 3 J J \z\ 2 D 2 a(\7a,VaJ dx dt + s 3 J J a\z\ 2 Aa\Va\ 2 dx dt 

fix(0,T) nx(0,T) 

„2 r r 

„2|„|2| a _|2 



a A \z\ z \Aa\ l dx. (2.23) 

ax(o,r) 

We finally estimate the terms involving the time derivatives in a: 

h = -\ jj M 2 (sd t V + 2s 2 d t (jVCT| 2 )) dx dt + s 2 JJ a\z\ 2 Aad t a dx dt. (2.24) 

nx(o,T) six(o.r) 

The next step is to give convenient estimations for the terms defined above. In order to do 
that several lemmas are proved. Making the notations ft ro — ft n 5(0, ro), O = f2 \ (ujq U ^ ro ): 
O = O \ fi ro we successively obtain 

Lemma 2.4. It holds that Im > 0, for any A > 0. 

Lemma 2.5. There exists Xq such that for any A > Ao and any s > then 
r r r „ - 2 

h > C 2 s 



f2x(0,T) 



|x| 2 ^|Vz| 2 + ^ 
1 11 \x\f 



dxdt + C 7 sA 2 // 6( l -^-) X 0\Vz\ 2 dxdt+ 



Ox(0,T) 



C*s\ 



•NUA-2 — 



— ) |Vz| 2 dx dt - B\s 
ro 

n ro x(0,T) Ox(0,T) 



C 6 sA 2 



\z\ 2 dx dt 



'^) A 0|Vz| 2 dxdt, (2.25) 

?*0 



w x(0,T) 

where C2,C7, C§,C§ and B\ are constants uniform in s and A, respectively uniform in s. 
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Lemma 2.6. There exists Ao such that for any A > Ao there exists so(A) such that for any s > sq 
it holds 



Inl > 



3 \x\ 2 \z\ 2 dx dt + C 15 s 3 A 4 / / 3 (l) 3 V 3 |x| 3A |z| 2 dx dt 

O ro x(0,T) Ox(0,T) 



C 1B s 3 A 4 



33/ 1 \3A -3i |3A I _ |2 



'r ' 



x\ 6 *\z\ z dx dt (2.26) 



w o x(0,T) 



for some constants C15, Ci6 uniform in s and A. 



Taking into account the negative terms in the expression of /; that we want to get rid of, we 
define 

I r =h- B x s ( f e\z\ 2 dx dt. 



(2.27) 



fix(0,T) 



Then 



Lemma 2.7. There exists Ao such that for any A > Ao there exists so(A) such that for any s > sq 
we have 



K r > 



W X \z\ 2 



Ox(0,T) 



dxdt+^ J J 6» 3 |x| 2 |z| 2 dxdt. 

Qt 

(2.28) 



Carleman inequality in the variable z 

From lemmas above we obtain 

Theorem 2.8. There exists a positive constant K and Xq such that for A > Ao there exists sq(X) 
such that for any s > sq we have 



sA 2 // 0(^) A 0|Vz| 2 dxdt + s 



TO 



H 2 -^|Vz| 2 



Ox(0,T) 



nx(o,T) 



\x\-r 



dx dt 



f!x(0,T) 

< K(sX 



3 3 A 4 [[ 3 (M) 3A 3| Z |2 dxdt 



Ox(0,T) 

e(M) A 0| Vz | 2 dxdt + s 3 A 4 
r 



6» 3 |:r| 3A c/> 3 |z| 2 dx dt). (2.29) 



w o x(0,T) w o x(0,T) 

Undoing the variables in (4.2) we obtain the conclusion of Theorem 2.2. 



3 Proofs of technical lemmas 



Lemma 3.1. Let us consider the function 



'0 



The Hessian matrix 



D 2 t^ 



3% 



dxidxj 



i,j=X,...N 
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is positive definite uniformly in \x\ < r$. for any A > 2. More precisely, for any A > 2 and any 
x G l w urai/i |x| < ro we Ziawe 



z \ro / 



Proof of Lemma 3.1. Firstly, we write D 2 t$(£,£) = <p(l/ro) x S^ where 

S,p : = A|C| 2 |x| A - 2 + A(A - 2)\x ■ £| 2 M A - 4 + 2A 2 (x • ■ W0M A ~ 2 + 
+ A|x| A ^V(e,e)+A 2 |x| A |VV;-C| 2 - 
Next we use to the inequality 



(3.1) 



(3.2) 



2\ 2 (x-Z)(Vtp-£)\x\ A ~ 2 <aA 2 |x-ei \x\ A_4 + — \x\ A \Vip ■ £\ 2 , V a > 

a 

Combining this with (3.2) we obtain 



> A|^| 2 |x| A - 2 + (A 2 - 2A - aA 2 )|x ■ C| 2 |x| A - 4 + X\x\ x D 2 ^, £) + (A 2 - ^)M A |W> ■ £| 2 . 

Next we choose a > such that A 2 — 2A — aA 2 = 0, that is a = (A — 2) /A. For this value of a we 
remark that 



o\2 

> A|£| 2 M A ~ 2 + X\x\ x D 2 ^,0 - -—I^IWIW, Vz,V£ 



A- 2' 



Then we deduce that 
A 



and finally 



> jM A - a |f| 9 + AM A -W(i - ^N 2 I^IL - \x\ 2 \d 2 mI 



s*>Z\x\*-*\tf 



(3.3) 



(3.4) 



holds true for A > 2 and for any |x| < ro. This yields the proof of Lemma 3.1. 

Proof of Lemma 2.4- It suffices to prove that Vcr • n > in dfl x (0, T). Firstly we have 

Vo- = 6{t) \-2xip- |x| 2 VV' - (— ) A (Ada?| A_2 + A|x| A V?/>)</> 
L V ' 

Dcpt to the first two conditions in (2.7) we have 

V^-n = -\Vij)\, WedfL 

In consequence we get 

Vct • fi = 6(t) [ - 2a?- fi+ |s| 2 |V^| + A(— )V|d A - 2 (b| 2 |VV>| - s • n) 
L v r y 



□ 



(3.5) 



Due to (2.8) we have 



Vcr . n > 0(f) |x| 2 (|VV| - 2C n ) + - C n ) 

L V 

which is positive since "0 satisfies the third condition in (2.7). This completes the proof of Lemma 
2.4. □ 
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Proof of Lemma 2.5. The computations here require a more careful analysis. First of all, we make 
the notations 

^2 = -6(t)T x 2, where t x i = \x\ 2 ip, 

respectively 

M a 

o~6 = -d(t)rs, where ta, = ( — ) <f>. 

and 

T = T x 2 + Ttf,. 

Next we write some formulas for t x 2 and that we are going to use in our computations. More 
precisely, for all x £ WL N and any i, j = 1, . . . ,N we have 

d Xz r x 2 = 2xiip + \x\ 2 d Xz ^ 
d 2 t x 2 = 2S. u i> + 2 Xl d X] i> + 2xjd Xi ip + \x\ 2 d 2 Xt V, (3-6) 



respectively 



d x ^ = (l) A (A^|x| A - 2 + X\x\ x d x ^)^ 
3U, T * = (-) A U% + A(A - 2)x l x \x\ x ~ i + X 2 x 3 d x M x - 2 + \ 2 Xi d Xj <ip\x\ x - 2 + 



I A— 4 , « „/,l„,|A-2 , \2„ ^ „/,l„,|A-2 

I AUj, T /\\/\ — LyjjiJGjXJu 

^0 

|Aq „/. i \2| lA; 



(3.7) 

Computations for 

Next we split in two parts as := 1} + if where 

1} := -2s J j D 2 a(Vz, Vz) dx dt - s JJ Aaa\Vz\ 2 dx dt + 2/is J j j^x • Vcr dx dt, 

f!x(0,T) Ox(0,T) Ox(0,T) 

J ' 2 := I // I z | 2a2 ' j ( 1 + a ) dxdt + s k| 2 Va • VAcr dx dt+ 

fix(0,T) Qx(0,T) 



,zpAcr Aa dx dt + (is JJ J^ Aa a dx dt ( 3 - 8 ) 

fix(0,T) £~2x(0,T) 

Moreover, we split as l\ := I^ x2 + l}^ where 
lj x2 = -2s II D 2 c7 x 2(Vz, Vz) dx dt - s 1 1 Aa x 2a\Vz\ 2 dx dt + 2fis 1 1 j^x ■ Vo^a dx dt 

(3.9) 



\x[ 

Ox(0,T) Qx(0,T) Qx(0,T) 



7^ := -2s JJ D 2 a 4> (Vz, Vz) dx dt - s JJ Aer a|Vz| 2 dx dt + 2fis JJ j^L x . \7o> dx dt. 

f2x(0,T) fix(0,T) ax(OJ) 

(3.10) 

Estimates for 2 : 
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Using the relations (3.6) we have 



Ox(0,T) 



|Vz| 2 -^j-j^JV dx dt + 8s J J 6(x ■ Wz)(Wip ■ Vz) dx dt 

Qx(0,T) 
|2 



-Is J J 6»|x| 2 L» 2 V(Vz,Vz) dxdt-2/is J J 6>j-^L(V0-x) dx dt-s ^ Acr^alVzl 2 dx dt. 

(3.11) 



nx(o,T) 



fix(0,T) 



f!x(0,T) 



Next, we estimate the first term in 1} 2 applying the Hardy inequality (1.5). We proceed as 
follows Firstly, fixing the time and integrating in the space variable we get the identities 



zVz • VV> dx = -| / z 2 A?/> dx, 



\ 2 -^zVz ■ W dx = — / |x| 2 - 7 A^ 2 dx- 



(2-7) 



(3.12) 



(x- V^)|x| _7 z 2 dx (3.13) 



Secondly, we apply the Hardy inequality (1.5) setting u := z\p§. Then, according to (3.12)-(3.13) 
and integrating in time we have 



Oipz 2 dx dt 



6 



|Vz| 2 - n— if; dx dt 



f!x(0,T) 



nx(o,T) 
1 

+ 4 

f!x(0,T) 



^^dxdt-i 

■0 2 



>C 2 



!ix(0,T) 

xr 7 ivzi 2 +^-i^dx dt+^ 



dx dt 



nx(o,T) 
2 



M 7 

|x| 2 ~ 7 A0z 2 dx dt - C* 2 



z 2 dx dt 



fix(0,T) 

(2-7) 



(x ■ vvoixpz 2 dx dt - 



nx(o,x) 



f!x(0,T) 



Since r satisfies (2.5) we have C 2 ip/{2\x[<) > C 2 (2 - 7 )|L>V|ooM 1 " 7 /2 for |x| < r . Then we 
easily obtain 



|Vz| 2 - M7 



V> dx dt > 



C 2 



x| 2 " 7 |Vz| 



ip dx dt 



fix(0,T) 



Ox(0,T) 



4 2' 

+ ^(2- 7 )r^ 7 |^U) JJ 6z 2 dxdt. (3.14) 

Qx(0,T) 

Denoting C 3 := dlVU + (AC/4 + |£ 2 0U/2 + Ga^pVloo^ + C 2 (2 - 7 )r^ 7 |^|oo/2 
and C 2 := C2/2, where i?n = su P xS n M> we 



\Wz\ 



M |xl 2 



•0 dx dt > C 2 



12-7 



|Vz| 



fix(0,T) 



fix(0,T) 



M 7 



ip dx dt 



6»|z| 2 dxdt. (3.15) 



Qx(0,T) 
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From (3.15) we obtain 



z 2 



x| 2 - 7 |Vz| 2 + f— dxdt-4sC 3 // 6>|z| 2 dxdt 
FpJ JJ 

Ox(0,T) Ox(0,T) 



-8s|DVU yy 6»|x||Vz| 2 dx dt - 2s|L»V|oo J J (9|:z;| 2 |Vz| 2 dx dt 

f!x(0,T) fix(O.T) 

-2s|/i||£)V|oo ^ ^^j-dxdt-s ^ Acr^alVzpdxdt. (3.16) 

fJx(0,T) Qx(0,T) 

Since r satisfies (2.5) we observe that 2C' 2 \x\ 2 ^ > SliPT/'loola;! + 2|Z) 2 -0| OO |x| 2 and 2C 2 /\x\~< > 
2\n\\Dip\/\x\, for all \x\ < r . Hence, from above we easily obtain 

ll x 2>2sC' 2 JJ 9 |x| 2 -t|Vz| 2 + dxdt -C 4 s yy 6»|z| 2 dxdt 

f2x(0,T) Ox(0,T) 



-s J J Aa x 2a\Vz\ 2 dx dt (3.17) 
nx(o,T) 

where C 4 = 4C 3 + 2|/x||£>V|ooAo. 
Estimates for : 

In order to get rid of the gradient term with negative sign in (3.17) we have to estimate from 
below the quantity T ■= I] a — s JJ Aa x 2a\\7 z\ 2 dx dt that is 

nx(o,T) 

T = -2s ([ D 2 a 4> (Vz 1 Vz) dx dt - s U Aaa\Vz\ 2 dx dt + 2^is 1 1 ^\-x ■ Vcr dx dt. 



Ox(0,T) Ox(0,T) Qx(0,T) 

To do that, according to Lemma 3.1, formulas (3.6)-(3.7) and (2.5)-(2.9), we firstly observe the 
pointwise inequalities 



2L> 2 Trf,(Vz,Vz) + AralVzl 2 > C 5 A(^) A 2 (/.|Vz| 2 , Vx e fi ro (3.18) 

|2D 2 r (Vz,Vz) + Ara|Vz| 2 | < C 6 A 2 (^) A </>|Vz| 2 , V x G oj . (3.19) 

Since (2.2) is accomplish we get 

2L» 2 r (Vz,Vz) + Ara|Vz| 2 > C 7 A(^) A 0|Vz| 2 , VxeO. (3.20) 

for some positive constants C5, Cq, C7 depending on tp. On the other hand, it holds that 

2 J^p± < C 8 X(^r% Vx€fl. (3.21) 
F r o 
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where Cg = Cs(ip,n) > 0. Therefore, it follows from above that 



T > C 7 sA 2 



<j)\Vz\ 2 dx dt- 



Ox(0,T) 

+ C 5 sX 



'M) A - 2 | V z| 2 dxdt-C 8S A 
' r 



) x -*<t>\z\ 2 dx dt 



n r x(o,t) 



Ox(0,T) 



C fi sA 2 



'M) A 0| V z| 2 dxdt (3.22) 
r 



w o x(0,T) 



Summing the terms in (3.17)-(3.22) we get 



I\ > C 2 s 



\x\ 2 ^\Vz\ 2 



\x\-y 



dx dt + C 7 sA 2 



r 



0|Vz| 2 dx dt- 



Ox(0,T) 



f!x(0,T) 

+ C 5 sA 

n ro x(o,T) ax(o,r) 

- C 4 s [ [ 6\z\ 2 dx dt - C 6 sA 2 



,^,A-2 2 dxdt _ CgsA 

" r 



^) A " 4 ^|z| 2 dxdt 

0(M) A 0| V z| 2 dx dt (3.23) 



nx(o,T) 



w x(0,T) 



Estimates for if . 

Making use of the support of a located far from origin we note that 

Ar 



|A 2 r*|, | At* |, |VAr |, 



< A\, VieSl, 



where A\ is a big enough constant depending on A. Then we get 

if > -a x s { t 9\z\ 2 dx dt 



Next we conclude 



h > C 2 s 



Ox(0,T) 



dx dt + C 7 sA 2 



^M) A 0| V z| 2 dx dt+ 
ro 



nx(o,T) 



Ox(0,T) 



+ C 5 sA 



■N\A-2 lr . 



|Vz| 2 dx dt - Sas 



I 2 dx dt 



n ro x(o,t) 



nx(o,T) 



C 6 sA 2 



(3.24) 



(3.25) 



u x(0,T) 



where B A = C 4 + A x + C 8 Asup xen {(|x|/r ) A </>}. 



^) A 0|Vz| 2 dxdt, (3.26) 
r 



□ 



Proof of Lemma 2.6. We split I n i = I n l,i + Inl,2, where l n i,\ are the integrals in I n i restricted to 
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fl ro and /„; 2 are the terms in /„/ restricted to O. We put a = —6t. Then I n i could be written as 



Inl = 2S 3 



3 \z\ 2 D 2 t(Vt, Vt) dx dt-s 3 J J 6» 3 |z| 2 aAr|Vr| 2 dx dt 

Qt Qt 



- — J J 9 2 a 2 \z\ 2 \Ar\ 2 dx dt. (3.27) 

Qt 

The term I n i is even more difficult to deal with because of technical computations. First of all we 
have 

8 Xi T = 2x t i> + \x\ 2 d Xi ip + (— ) A (A^|a;| A - 2 + X\x\ x d Xi i(i)<t>, 
d 2 tXj T = 26^ + 2 Xl d Xj yj + 2 Xj d Xi i) + \x\ 2 dl. x .ip + (^-) A [a^ + A(A - 2)x i x j \x\ x - 4 + 



+ \ 2 x 3 d x MA + x 2 x,d x ^\x\ A - A + \\x\*d XiX .i> + \ 2 \x\*d Xi i)d Xj i> 

Using the expressions in (3.28) we obtain several useful formulas: 

\x ■ Vr| 2 = M 2 |Vr| 2 + (|V</> ■ x\ 2 - \x\ 2 \^\ 2 ) (\x\ 2 + A(1) A ^xrf \ 
(x ■ Vr)(W ■ Vr) = |Vr| 2 W ■ x + (|x| 2 |Vi/;| 2 - \x ■ V^| 2 ) (l + A(— ) A 0|cc| A ^ 2 ) : 
x (2VN 2 + A(1) A 0N A ), 
|W • Vr| 2 = |VV^| 2 |Vr| 2 + (\x • Wf - M 2 |W| 2 ) (|x| 2 + A(^) A 0|x| A ) ' 
Using the expressions in (3.29) we reach to the identity 



(3.28) 



(3.29) 



2D 2 t(Vt, Vt) - aAr|Vr| 2 = 

= 2^(2 - aA0|Vr| 2 + 4(2 - a)|Vr| 2 VV> ■ x + 2\x\ 2 D 2 ^(\/t, Vt) + 

+ 8(N 2 |V^| 2 - |VV • .*| 2 )(l + A(i) A ^N A - 2 ) fa\x\ 2 + A(l) A 0| a; | A )- 
- a\x\ 2 ^\Vr\ 2 + {^r) X (f>{ [((2 - a)A 2 - A(2 + aN - 2a)) M A ~ 2 + 

+ 2A 2 (2 - a)|x| A " 2 VV -x + (2- a)A 2 |x| A+2 |V-0| 2 j |Vr| 2 + 
+ (2 - o)\\x\ x D 2 ^{Vt 1 Vr) + (|x| 2 |V?/>| 2 - \Vtp ■ x\ 2 )x 



1 \ 2A 

'TV 



4A 3 (-)'> 2 |.t| 3A - 4 + 8A 2 (-) A 0|x| 2A - 

r ' 



+ 8A 2 V(1- VON* -2A(A-2)|a;| A }. 



Computations for I n i,i- 

Taking into account the properties of ip and a and Lemma 3.1 we can obtain 



2L> 2 t(Vt, Vt) - aAr|Vr| 2 > |Vr| 2 + A 2 



1 nA 



-r ' 



x| A - 2 |Vr| 2 , Vi€f] r 



(3.30) 



On the other hand we have 



C 9 \x\ 2 > |Vr| 2 > \x\ 2 , \/xeQ ra , 
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for some constant Cg = Cs(ip) > 0, and therefore 



2 j D 2 t(Vt, Vr) - aAr|Vr| 2 a 2 > |x| 2 , VieSJ 



»T) - 



a 2 |Vr| 2 < C 9 |x| 2 , Vxen ro . (3.31) 

In particular we obtain 

Inl,i>Y jj ^M'M'dxdt. (3.32) 

n ro x(o,T) 

for ,s large enough. 
Computations for I n i t 2- 

Again, according to the properties of ip and a we observe that 

2D 2 t(Vt,Vt) -aAr|Vr| 2 > Ci A 2 ( — ) X <p\x\ x |Vr| 2 

>C u A 4 (-) 3A |x| 3 V 3 , Vx-eO, 

a 2 \Ar\ 2 < Ci 2 A 4 (-) 2A 2 |a;| 2A , VxeJ] (3.33) 
|2D 2 r(Vr,VT) -aAr|VT| 2 | < C 13 A 4 (-) A 0|x| A |Vr| 2 

< Ci 4 A 4 (-) 3A |x| 3A 3 , Vz€w (3.34) 
for some uniform constants Ciq, Cn, C12, C13, C14 depending only on ip. Then we obtain 



jj e 3 (i) 3 VN 3A k| 2 dxdt-C7 14 s 3 A 4 |y 3 (l) 3A ,/) 3 |x| 3A |z| 2 dxdt- 



^ni,2 > CioS 3 A 4 

Ox(0,T) w x(0,T) 



Ci2.s 2 A 4 /Y 2 (l) 2A |x| 2A ^ 2 |^| 2 dxdt. (3.35) 

Ox(0,T) 



By summing the terms in (3.35)-(3.32) we obtain 



/„;>y yy ^i^Hzi 2 dx dt + c 15S 3 a 4 yy e 3 ( ; -) j > 3 N 3A izi 2 dxdt 



n ro x(0,T) Ox(0,T) 



-C 16 s 3 A 4 yy 3 (l) 3A 3 |x| 3A |z| 2 dxdt (3.36) 
where C15 = Cio/2 and Ci6 = C12 + Ci4- 

□ 

Proof of Lemma 2. 7. According to the expression of 9 we obtain 

\e'\ < ce 1+1 ' k , \e"\ < ce 1+2 ' k , 

for some positive constant C . On the other hand according to the definition of a we get 

|A<t| < D x e, v x e 0, 

|3eo-| < D\6' , V z e n, 
5 t (|V<7| 2 ) < D x 66'\x\ 2 , V ief! r( , 

5 t (|VfT| 2 ) < D x ee'\x\ 2X , v ied, (3.37) 
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for a big enough constant D\ > depending on A. Since a is supported far from origin we can 
write 



a\z\ 2 A<rd t <T 



AD 



dx dt < — =±s 



A „2 



e 2+1 i k \x\ 2 \z\ 2 dxdt 



£2x(0,T) 



a ro x(o,r) 



£ 2 S 2 I f e 2+1 ' k \z\ 2 dx dt (3.38) 



Ox(0,T) 



From (3.37) we also obtain 



||z| 2 9 t (|Va| 2 )| dxdt<C A s 2 jj e 2+l / k \x\ 2 \z\ 2 dxdt 



nx(o,T) 



Now, we put 



n,. x(o,t) 



+ D\s 2 JJ 6 2+1 l k \x\ 2X \z\ 2 dxdt. (3.39) 

6x(0,T) 



i?:=-- |z| 2 9 tt( r dx dt - B x s 



\z\ 2 dx dt, 



fix(O.T) 

where B\ is elected as in Lemma 2.5. 
Then, there exists E\ > such that 



nx(o,r) 



1 2 df* (j dx dt 



< 



E x s JJ 9 l+2 / k \z\ 2 dxdt 



Qx(o,T) 



Ox(0,T) 



B A s y ^ 6»|z| 2 dx dt < E x s JJ 9\z\ 2 dx dt. 



(3.40) 



fix(0,T) 

Summing these bounds we obtain 



Ox(0,T) 



\R\ < 2sE x JJ 6 1+2 ' k \z\ 2 dx. 

f!x(0,T) 



(3.41) 



In the sequel, we write 



J j 6 1+2 ' k \z\ 2 dx dt = JJ (>i+2A-W \ x \ih' (I^Vff'j a .|-y/«'|«|a/9') dx dt. 



Qt 



Now we take 



2+7 , 7+2 

g = , q - 



2 

Note that 1/q + 1/g' = 1 and applying the Young inequality we obtain 

V+ 2 / fc H 2 dxdt</3« ^ ^ 1+2 / fc - 1 /9') < ?|x| 2 |z| 2 da;dt+ 7 i 7 ^ ^ dx dt 

Qt Qt Qt 

= fJJ 3 M 2 M 2 dx dt + ip yy^ dxdt, 

Qt Qt 



(3.42) 
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provided k = 1 + 2/7. Therefore, from (3.41) and (3.42) we have 

\R\ < C{\)a(p f f 3 M 2 |z| 2 dx dt + -L // (?M1 dx dt) . (3.43) 



N_ 2 

Qt Qt 

Consequently, from (3.38)-(3.43) follows 

\I T \<F x (s 2 JJ e 2 \x\ 2 \z\ 2 dxdt + spo JJ 6 3 \x\ 2 \z\ 2 dxdt 

£V X(0,T) Qt 

2 dxdt + s 2 // 6» 3 |a;| 3A |z| 2 dx dt ) , (3.44) 



P q ' J J M 7 

nx(o,T) <5x(o,T) 

for a new constant F\ > 0. 

Take j3 such that F\/ (3 q = C2/2. Then there exists some s (X) such that for s > s (X) we 
finish the proof of Lemma 2.7. □ 



4 Sketch of the proof of Theorem 1.2 



Let us also assume w C O is a non-empty control region such that Q n B\ (0) is included in fl and 
£?i(0) Ho) is empty. This can always be done by a scaling argument. Moreover, let us also consider 
a non-empty open subset loq CC u>. 

Again, the proof is based on Carleman estimates as shown in Theorem 1.1. Before sketching the 
proof it is worthy to make some remarks about the technical difficulties that appear for domains 
Q not satisfying the property (P) (as in the hypothesis of Theorem 1.2) when choosing the weight 

If it were possible to construct a weight ip as in Theorem 1.1, for a domain which does not 
satisfy the property (P), then the result of Theorem 1.1 could be generalized to any geometrical 
configuration of f2 in the sense that we could control our system with a distributed control in any 
open subset wCfino matter what the geometry of Q is. Unfortunately, for a domain f2 as in the 
hypothesis of Theorem 1.2 we have not found any example of such ip as in Subsection 2.1, since 
the distance function d = d(x, dCl) mentioned in Remark 2.1 violates the condition (2.9) in this 
geometrical configuration. 

Despite of this, for the proof of Theorem 1.2 (ft not satisfying the property (P)) we choose the 
weight ip to satisfy the same constraints as in Subsection 2.1 less the condition (2.9). In view of 
that, the result of Theorem 1.2 is weaker than the result provided by Theorem 1.1. Here we skip 
the details of the proof of Theorem 1.2 since it is straightforward from the Carleman estimates 
developed for Theorem 1.1, pointing out just the conclusions. 

Indeed, when doing Carleman estimates, because the loss of condition (2.9), we cannot capture 
the energy in the annulus region 

J := {x e | r < |x| < rj, (4.1) 

where r$ stands for the constant in (2.5) (ro is small enough since it is proportional with 1/|V?/>|) 
and n is any fixed positive constant such that r± < 1. For that reason, we succeed to control 
our system forcing the control to act precisely inuU lo' , where ui' is defined in (4.1). If lo is a 
non-empty open subset w CC w' we obtain the following Carleman estimate: 

Theorem 4.1. There exists a positive constant K\ and Xq such that for X > Xq there exists sq(A) 
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such that for any s > sq we have 



e- Zsa \Vw\ A dx dt + s 



-2sa 



x| 2 - 7 |Vu.| 2 + 



Ox(0,T) 



fix(0,T) 



dx dt 



33 e- 2s,7 |x| 2 | W | 2 dx dt + ^A 4 



3\4 f f a 3-2sa( M \ 3 > j.3|„„|2 



wr dx dt 



Ox(0,T) 



< if sX 2 



( sA 2 yy 0(M) A 0e -2-| Vw | 



Vtul dx dt + s^A 



©x(0,T) 



3\4 



? 3/Nn3A ,3 e -3«n |2 dx dt 



r 



w Uii; x(0,T) 



w o UWoX(0,T) 



(4.2) 



In particular, according to the Cacciopoli's inequality, Theorem 4.1 yields the proof of Theorem 

1.2. 
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